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(— i , The flavour violating neutrinoless muon decays in the presence of nuclei, are discussed. 

■^j" \ We focus on the theoretical aspects of u~(A, Z) — > e~(A, Z)* (muon to electron conversion) , 

one of the most prominent flavour changing reactions, emphasizing its connection with 
the physics beyond the standard model. This process offers the most severe limits of the 

\Q \ lepton flavour violation. By using the nuclear transition matrix elements calculated with 

several methods, and the recent experimental data of the branching ratio Rn e - , we determine 
limits for the flavour changing parameters entering the elementary sector part of R^ e - ■ 
These results are discussed in view of the ongoing experiment at PSI and the designed at 
Brookhaven, which are expected to push down by some orders of magnitude the experimental 
sensitivity the next few years with the hope to see "new physics". 
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1 Introduction 

In the framework of the standard model of the electroweak interactions, if one associates with 
each lepton generation a lepton flavour quantum number as 

(for antiparticles L e , L M , L T take values with opposite sign), the quantum numbers L e , L^, L T 
are separately conserved. In addition, the total lepton number 

L = L e + L^ + L T (2) 

appears also to be conserved 0-0 • Even though these conservation laws are uncontradicted with 
experiment || ||, the standard model is believed to be an effective low energy approximation 
of a more fundamental theory which is unbroken at high energies. In this spirit, many theories 
beyond the standard model have been developed during the last decades (grand unified models, 
supersymmetric theories etc.) pJ-fOJ. in which the existence of flavour and lepton changing 



effective currents is predicted. In this way, a plethora of exotic processes, non-conserving the 
flavour, L e , L M , L T , and/or lepton, L, quantum numbers, are possible. For some examples and 
the existing limits see Table 1 @, |?|-|T5|. 



In the present work we are going to discuss the category of neutrinoless exotic muon decays 
which take place in a muonic atom, i.e. in the case when a muon is bound in the field of a 
nucleus. In the context of the aforementioned extensions of the standard model this bound muon 



can give two yet non-observed processes; the conversion of \i b into an electron [[R|-|32 



u b + (A, Z) — ► e + (A, Z)* (fj, b — ► e conversion) (3) 

which violates the lepton flavour quantum numbers L e and L M , and the conversion of the muon 
fi b into a positron 



a b + (A, Z) — > e + + (A, Z — 2) (fi b — > e + conversion) (4) 

which violates the conservation of the total lepton L as well as the flavour L e and L^ quantum 
numbers. 

The important role played by the muon as a test particle in the investigation of the electro- 



magnetic and weak interactions with hadronic particles, has been recognized decades ago [33 



The bound in the innermost Is atomic orbital muon \i h has especially prominent role, since it 



gives rise to very important allowed decay modes |34| , |35|[ . it can disappear either by decaying 
into an electron and two neutrinos, 

ji b — > e~ + v e + v^ (muon decay in orbit), (5) 



or by being captured by the nucleus |34|, [35 



fi b + (A, Z) — > v^ + (A, Z — 1) (ordinary muon capture) (6) 

The latter two processes conserve both lepton and flavour quantum numbers and have been well 
studied theoretically as well as experimentally |34|, |35 . 



The exotic channels of Eqs. (^|) and (|4]), are the subject of some of the most important recent 
experiments ]5|, || searching extensively for events involving muon number non-conservation. Up 
to now, these experiments have only put upper limits on the branching ratios R^ e - and -R Me +, 
i.e. on the rates of processes (^) and (||) divided by the total rate of the ordinary muon capture 
reaction. For a comparison of these upper limits with those obtained from other lepton flavour 
violating reactions see Table 1. Obviously, the observation of any of these processes will signal 
"new physics" beyond the standard model and restrict the various models predicting such exotic 
reactions. 



1.1 Basic features of the exotic muon decays in nuclei 

The neutrinoless rare decays, ^b — > e~ and fib — > e + conversions in nuclei, start from the same 
initial state (a muon at rest in a muonic atom). Experimentally, they are both simultaneously 
studied. Theoretically, both of them can occur by assuming mixing of intermediate (massive) 
neutrinos, but (/i~, e~) cannot distinguish between Dirac and Majorana neutrinos while (fi~ , e + ) 
can proceed only if the neutrinos are Majorana particles. These processes can also be predicted 
by supersymmetry via the mixing of the intermediate supersymmetric particles, e.g. s-leptons 
(see sect. 3). 

There are, however, big differences between the exotic decays (fi~,e~) and (/i~,e + ) origi- 
nating from the individual characteristics of their transition operators. Thus, due to its charge 
conserving operator, the Hb — > e~ conversion can take place in one nucleon; it can proceed via 
the coherent mode (when the participating nucleus remains in its ground state and all nucle- 
ons participate) JO]], the only measurable (fi~,e~) conversion channel, which is an important 
experimental feature. 

Another experimental advantage of the (/z~, e~) we should mention is the fact that, a detection 
of only one particle is sufficient and there is no need for coincidence measurement M. If the 
electron energy is equal to E e « m^c 2 — e^ (coherent mode) the process is free from the reaction 
induced background which involves, in addition to the muon decay in orbit mentioned above, the 
following reactions: 

i) Muon decay in flight. 

IT -> Uf, + v e + e~ (7) 

ii) Radiative muon capture followed by the creation of a fully asymmetric e~ , e + pair at zero 
energy, i.e. 



ft + (A,Z)^(A,Z-1) + u^+ 7 (8) 

1 -> e + + e" 

On the other hand, the muon to positron conversion (/i~,e + ), is more complex since it is 
a two-nucleon process; it is similar to but more complicated than the neutrinoless double beta 
decay 

(A, Z) -* (A, Z±2) + e T e T (Ovpp - decay) (9) 



which is a lepton and flavour number violating process extensively discussed in the present 
workshop (see e.g. pUf). 



Because of its simplicity and important characteristics, the (/i , e ) conversion is preferentially 
studied by many authors [f7|-|I3|, and in the present work we will devote to it a rather extensive 



discussion. We pay special attention to the construction of the effective Hamiltonian in the 
framework of several common extensions of the standard model. This is subsequently used 
to calculate the relevant transition rates. Using the results of the branching ratio R^ e - and 
exploiting the existing experimental data we put limits on the flavour violating parameters (see 
sect. 6). 



2 The (/i , e ) conversion in Nuclei 

2.1 Brief historical review 

The first experimental searching for (/i~,e~) events has been done long ago by Steinberger and 



Wolfe [[16) at Columbia Cyclotron using Cu as a target. This experiment put an upper limit for 
the branching ratio -R Me - 

R,e~ = w r( ^" , (10) 

T{fj,~ — > capture) 



on the value -R Me - < 5 x 10 . Some years later Conversi et al. |17j in the synchrocyclotron 
at CERN using again Cu as target with two independent experiments found the limits R^ e - < 
5 x 10~ 5 and R^ e - < 5 x 10~ 6 . A Cu target was also used at TRIUMF by Bryman et al. (18 
in the TPC experiment. From the data of this experiment the limit -R Me - < 2.6 x 10~ 8 was 
extracted. 

We should note that, the use of a Cu target in the aforementioned experiments had been 
influenced by the theoretical estimates of Weinberg and Feinberg [0] who found that, the ^ — > e~ 
conversion branching ratio shows a maximum in the region of Cu (see also Ref . [^5| ) . 

Target different than Cu was first used at SIN (present PSI) by Badertscher et al. |19| in the 
SINDRUM spectrometer who used the isoscalar target S. The upper limit set on R^e- by this 
experiment was R^ e - < 1-6 x 10~ 8 . 



In another TPC experiment at TRIUMF (Bryman et al. |S{J) a Ti target was used and pushed 
the upper limit of R^ e - at first to the value -R Me - < 1.6 x 10~ n and some years later [21] to the 



value R^- < 4.6 x 10~ 12 . 

The best upper limit on i? Me - set up to the present has been obtained by using Ti target at 
PSI which yielded the value R^ e - < 7.0 x 10~ 13 . Recently, by using a heavy 208 Pfr target in 
the SINDRUM II experiment the limit -R Me - < 4.6 x 10~ n was determined 0. This limit was 
an improvement by an order of magnitude over the previous limit R^e- < 4.9 x 10~ 10 obtained 
from preliminary results on Pb at TRIUMF (Ahmad et al. f21~H ). 



Ongoing experiments at PSI as well as planned experiments at Brookhaven, discussed by 
Walter || in this conference, are expected to increase the experimental sensitivity by two to 
three orders of magnitude with the hope to see some events of the reaction @ in the near future 
or to drastically push down the best upper limit of the branching ratio R^- . 



From a theoretical point of view, the basic background for the reaction ^ —>■ e~ has been 
formulated by Weinberg and Feinberg [[!]] who assumed that this process is mediated by a virtual 
photon (photonic mechanism). Non-photonic contributions were included later on in the post- 
gauge-theory era |H|-IT3|. Afterwards, a lot of new mechanisms which lead to the \i — e conversion 



have been proposed |30[. The mechanisms which involve intermediate neutrino mixing or mixing 
of intermediate s-fermions are extensively discussed in the next section. 

From a nuclear physics point of view, it is important to know the transition matrix elements 
of the participating nucleus in processes (0) and (|]). For the reaction (^) this was recently done 
for a series of nuclear isotopes |31| including 48 Ti and 208 Pb which are of current experimental 



interest || p3| . Several models, like shell model, local density approximation, quasi-particle RPA 
etc., have been employed (see sect. 4). 



Effective Hamiltonian of (// , e ) in common extensions 
of SM 



Some Feynman diagrams, which contribute to fi^ — > e~ conversion in the lowest order of per- 
turbation theory (one loop level) are shown in Figs. 1 and 2. In Fig. 1 mixing of intermediate 
neutrinos and gauge bosons is involved and in Fig. 2 the mixing of neutralinos and s-leptons in 
a supersymmetric model is considered. In the context of these models we construct below the 
effective Hamiltonian for the (/z~, e~) conversion by first writing down the hadronic and leptonic 
currents of these diagrams as follows. 



3.1 Hadronic currents at the nucleon level 

Usually one starts from the weak vector and axial vector quark currents. By assuming a specific 
nucleon model, one transforms these hadronic currents to the nucleon level. For the diagrams 
shown in Figs. 1 and 2, one gets the following results. 

i) Photonic diagrams (Figs. 1(a) and 2(a)): In this case the hadronic vertex is the usual 
electromagnetic coupling 

J« = N plx N p = N lx l - (1 + r 3 ) N (photonic) (11) 

[N \ 
N is the nucleon isodoublet, N = I AT P ), with N p , N n the proton, neutron spinors. 

ii) Non-photonic diagrams: There are many such diagrams contributing to the /i^~ — > e~ 
conversion. The most obvious are: a) those mediated by Z-particle exchange (see Figs. l(a),(b), 
2(a),(b)) and b) the box diagrams involving either massive neutrinos and W-bosons, as in Fig. 
1(c), or s-leptons and neutralinos, as in Fig. 2(c). The non-photonic hadronic current which 
includes all the above cases can be compactly written as 



Jx ] = Nlx\ [(3 + f v (3r 3 ) - (f v p" + f A (3'r 3 ) l5 \ N (non - photonic) (12) 



where the parameters /?, /?', /3", for the models assumed take the values 
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(13) 



(14) 



(15) 



The parameter j3 is defined as (3 = /3i//3o, with @ (fti) the isoscalar (isovector) couplings at the 
quark level and /y, Ja represent the vector, axial vector static nucleon form factors (/a//v = 
1.24). We should note that, in general, the parameters (3, j3 and j3 are functions of sin 2 8w 
For example, in the case of Z-exchange we have (3 = 3/2sin 2 8w = 6.90 (for some other cases see 
Ref. 10011). 



3.2 The leptonic currents 

The corresponding leptonic currents in the diagrams of Figs. 1 and 2 are given from the following 
expressions. 



i) photonic mechanism: 



\v Qv 



J(i) = M ^ e ) ^ Ml + lsfEi)io- v — + (Ieo + 7s/mo) lv 9 



j* ^^ 



m. 



u(Pt> 



(16) 



ii) Non-photonic mechanism: 



J(2) = u (Pe)l X -z (A + /275)«(Pm) 



(17) 



where q = p e — p^ is the momentum transfer with p e , p^ the lepton momenta. The parameters 
Ieo, Iei, fhio, /mi, /i, /2 are gauge model dependent form factors. In the models mentioned 
above they are interelated as follows. 



3.2.1 Models involving intermediate neutrinos and gauge bosons 



For purely left-handed theories we have 

IeO — ~ /MO, JEl — ~ /Ml, 

In models involving VT-bosons, it, furthermore, holds 



/go 

q 2 



/i = -/ 2 



mf 



(18) 



(19) 



3.2.2 Models involving intermediate neutralinos and s-leptons 

In this case, the photonic and non-photonic parameters are related as 

(47ra)/ M i = ~{^a)f El = ~f (20) 

(4ira)f E0 = -(47to)/ mo = ~f (21) 

h = -h = ^/V, A) = -\sin 2 e w (22) 

lb 6 



. 9 '/? '/ - 9 

1 ^ TOO ^ 



/ = « 2 -f r/, fj = ^^P (23) 



m? are the average of the square of the s-fermion masses and n depends on the details of the 
model. For Z-exchange contribution in SUSY models (see also Ref. |3"0|] ) we have 



„2 m L 



h = -h = &^f, (24) 



m 



£ depends on the specific model; it vanishes if one ignores the Higgsino components of the 
neutralinos (e.g. in the case of a pure photino); it also vanishes if the probabilities of finding 
the two Higgsinos in the neutralino are equal. In the model 2 of Ref. |¥T] it takes the value 



£ = 2.8 x 10 2 . In models in which m 2 is larger than m|, the Z-exchange is more favored. Due 
to the smallness of £, however, many authors consider Z' etc. ||12| . 



3.3 The (/i , e ) conversion effective Lagrangian at nucleon level 

The effective Lagrangian at nucleon level, which contains contributions from all types of diagrams 
mentioned above, is written in the form 

KA .. 47m ,-A 7(1) , C A x(2) f2 ,x 

M ~ ~^r^w J x + ^" j (2) j a W 

(Photonic) (Non — photonic) 
where a is the fine structure constant, and 

t — J F yFi " ' ^ — boson exchange ,„„■< 

1, s — lepton mixing 



3.4 Transition operators at the nuclear level 

In order to obtain the transition operators at the nuclear level, where the calculations of the 
transition rates are done, one usually assumes that the lower components of the nucleon spinors 
is negligible. This assumption gives the non-relativistic limit of hadronic currents. For a nucleus 
(A, Z) this leads to the following fi^ — > e~ conversion operators 



^o = 9v E ( 3 + /WH-) e " Q ' rj (Vector) (27) 

3=1 

which is of Fermi-type, and 



A 

n = ~~9aY, (fvP" + fAP'r 3j ) -j= e"*™ {Axial Vector) 



(2? 



which is of Gammow- Teller type (fy = 1.0, J'a = 1.24). The parameters gv and g^ take the 
values gy — 1/6, g^ = 0, for the photonic case, and gy — g^ — 1/2, for the non-photonic case. 

In the calculations of the nuclear part of the transition rate, the following multipole expansion 
operators, obtained from Eqs. (|27|) , (|28|) , enter. 



f£ 0)J = <5 ZJ ^£(3 + fvPr^jiikrjYUti) (29) 



S(i,0)J 

— <JU V -±'i 

f£ 1)J = v^ £(/v/3" + /a^tw^CAtO^^Ox^]' (30) 

i=l 

The matrix elements of these operators calculated in a given nuclear model provide the indi- 
vidual multipole contributions (see Refs. [ F29] , |5T] , [52]]). From kinematics one can show that the 
momentum transfer q can be approximated by 

q =q f = jrift - e b - (E f - E gs ) (31) 

where e& represent the muon binding energy and (Ef — E gs ) the excitation energy of the nucleus 

(A,Z). 



4 The (/i , e ) conversion branching ratio R^ e - 

In studying the \i~ — >■ e~ conversion process, the most interesting quantities both theoretically 
and experimentally are (i) the branching ratio R^ e - defined in Eq. ([K]) and (ii) the ratio of 
coherent to total rate 

V = j^f tt 32 

In general, it is impossible to separate the nuclear structure dependence of R^ e - from the 
elementary particle parameters. In the coherent mode, however, this is sometimes possible (see 
Ref. [j3U|] ). In the models mentioned in the previous section the dominant coherent channel, can 
be written as 

Rfie- =91 (33) 

where p is independent on nuclear physics. The function j(A, Z) contains all the nuclear infor- 
mation and is defined as 

EePe M 2 

7 '" ml G*Zf GP (A,Z) { ) 



(G 2 ~ 6) with fop the Goulard-Primakoff |34j] function describing the total p^ — > v^ rate and 
M 2 the nuclear transition matrix elements of the \x^ — > e~ conversion. Thus, the nuclear aspect 
of the p — e conversion process involves the matrix elements M 2 which in general are given by 

M 2 = f 2 | (/ | Q | *,/i> | 2 +3/ 2 I (/ I ft I 2,/i) | 2 (35) 

where |/) the final nuclear state populated during the process. 



The quantity r\ of Eq. fl3"2|), in the case of the approximation Eq. (^) can be written in terms 
of the matrix elements of Eq. (EB) as 



T coh (Li^e-) . M ( 



2 
coh 



rj = -p^: Pf « =-£* (36) 



r toi (/i -> e~) M r 



tot 



where M 2 oh , M 2 ot the coherent, total rate matrix elements, respectively. Such calculations of r) 
are discussed in sect. 6. 



4.1 Definitions of j(A, Z) and p in various models 

i) Neutrino mixing models: For the photonic diagrams, the function 7 of Eq. (|4j) takes the form 

EePe Z\F z (q 2 )\ 2 
lph ml G 2 f GP (A,Z) [6() 

For the photonic mechanism with left-handed currents only, the parameter p of Eq. (|33|) is given 
by 

1 a \2I/ati + /eo| + |/ei + /mo| 9a m e 2 

p = (47ra) (g^p = 6i^<^ +wl (38) 

For the non - photonic diagrams the function 7 is given by Eq. (133]) with M 2 given by 



M 2 



1 + 3-f v PN F N (q 2 



3 + f v P Z F z ( 



lph (39) 



(Fz(q 2 ), F/v(g 2 ) are the proton, neutron nuclear form factors) and the parameter p associated 
with intermediate neutrinos is given by 



|A)/l| 2 + |A)/ 2 | 2 9 2 x2,on m e , 12 ,„ n x 

The quantities ??„ and r\ jv in Eqs. (|38|) and (|40|) are the flavour violating parameters associated 
with intermediate light (z/j) or heavy (Nj) neutrinos. They depend on the gauge model (we have 
ignored the neutrino mass independent contribution arising from the non-unitarity of U^) as 

3 e 

,-E^ef (-««f + a) ,4 2) 



Ugj , U^j , with i = 1 for light neutrinos and i = 2 for heavy neutrinos, are the elements of the 
charged lepton current mixing matrix associated with the electron and the muon as 

n n 

^ = E^S + E^i ( 43 ) 






"» -E^S + EO^ (44) 



where z/ e , ^ are the weak neutrino eigenstates. The number n counts the generations assumed. 

ii) Neutralino and s-lepton mixing models: Under some plausible assumptions and neglecting 
the Z-exchange, the function 7 and the parameter p are given by 



7 = Cbr + oVi^) i» ( 45 ) 



•13 1 N F N \2 
^12 + 2 Z ^ 

1 a 4 



P == 288 (^P ^ (46) 



In the case of Z-exchange diagrams we have 

- 14.04 

Z F z 



7 = 0.053(1- 14.04^ %P)\ h (47) 



and the quantity p is now given by 

= R R x 1 0~ 4 

(Gprn 



1 

-1 a 1-12 



5-5 x l(T 4 77r -^ If/I' (4^ 



We should note that, in the models discussed above, the only variables of the elementary 
particle sector which enter the function 7 are the parameters (3 = (3i//3 , (3' and (3" . Once 
7 (A, Z) is known, e.g. by nuclear model calculations, from Eq. (^) one can extract information 
about the interesting parameter p from the experimental data (see sects. 5 and 6). 



5 Nuclear model calculations for the branching ratio R^ e - 

The main aim of the nuclear calculations of the pT — » e~ conversion is to study the nuclear 
structure dependence of its coherent, incoherent and total rates throughout the periodic table 
and find favorable nuclear systems to be used as targets in searching for lepton number non- 
conservation. Furthermore, by exploiting the nuclear part of the branching ratio in Eq. (^) we 
can put restrictions on the elementary sector part of R^, e - , i.e. the quantity p, which contains 
the lepton flavour violating parameters (neutrino masses, mixing angles etc.) of a specific gauge 
model. In this way, one can check the models predicting flavour non- conservation and improve 
them. 



5.1 Calculation of the muon-nucleus overlap integral 

As we have seen in sect. 4, by neglecting the effect of nuclear recoil, the nuclear dependence of 
the p — e conversion rate is, in general, included in the matrix elements of Eq. (|35|) . 



In the special case of the coherent process, i.e. ground state to ground state transitions, only 
the vector component of Eq. (KH) contributes and one obtains 



(/|fio|^> =g v (3 + f v f3)F(q 2 ) (49) 

where F(q 2 ) is the matrix element involving the ground state proton, p p (x), and neutron, p n (x), 
densities (normalized to Z and N respectively) as 

F(q 2 ) = W) + |^k(g 2 ) (50) 

<J + JvP 

with 

F p ,n{q 2 ) = [d 3 x p p , n (x) e~^ x $„(x) (51) 



Thus, the coherent rate matrix elements can be obtained from Eq. fl49| ) by firstly calculating 
F(q 2 ) from Eq. (|50"D for a given muon wave function and given proton, neutron nuclear density 
distributions needed in Eq. (|5"T|). 



The methods used up to now, according to how do they treat the muon-nucleus overlap 
integral appeared in the branching ratio, are classified in two categories: i) those using exact 
evaluation of the muon wave function, and ii) those using the factorization approximation for the 
muon wave function. In the next subsections we present a brief description of these methods. 



5.1.1 Exact Evaluation of the muon wave function 

For point like nuclei the muon wave function $^(r) is trivially obtained by solving the Schrodinger 
or Dirac equations for the Coulombic nuclear potential. By taking into consideration the effects 
of finite nuclear size and vacuum polarization one needs to solve numerically these equations. 



Such calculations were done in Refs. p8| , [42| . In Ref. |28[ the muon wave function was obtained by 
solving numerically the Schrodinger equation taking into account the aforementioned effects. In 
Ref. |42| the Dirac and Schrodinger equations are solved using modern neural network techniques 



and an analytic expression for $ M (r) is obtained by means of optimization methods. 

5.1.2 Factorization approximation for the muon wave function 

A reasonable assumption for the Is muon wave function of light and medium nuclei is to assume 
that it varies little inside such systems. In this case one can use for the matrix elements Eq. fl4"9|) 
the following approximation 

(/ | fio | i, fi) « M3 + fvP) ($!.)</ | fio | *) (52) 

where ($i«) is the mean value of the muon wave function given by 

($ ls ) 2 = a -*U- (53) 

7T Zj 

a is the fine structure constant and Z e ff represents the effective charge which sees the \i^ in a 
muonic atom. For light nuclei /i^ lives outside the nucleus, Z e ff w Z, but for heavy and very 
heavy nuclei the \i^ resides inside the nucleus and Z e ff <^. Z. 



Using this approximation, for the coherent rate many authors (see Ref. [0, [11], p5| ) calculated 
the quantity F(q) in eq. (^) from the equation 

|F( g 2 )| 2 « ^< % |ZF z (g 2 ) + |^M iVi^g 2 )! 2 (54) 

7T Z d + JvP 



The nuclear form factors Fz, Fn are either calculated by using various models, i.e. shell model |25 



43| , quasi-particle RPA |29| etc., or extracted from experimental data whenever possible jP|, |4"5 | 



The branching ratio R^ e - in this approximation takes the form of Eq. ( |55| ) with 7 given as 
explained in sect. 4. 



5.2 Evaluation of the partial \i — > e conversion rates 

The methods used up to now, according to how do they evaluate the partial rates, i.e. the 
transition rate of each individual channel, are classified in two categories: 1) those doing state-by 
state calculations of the partial rates 2) those employing an effective evaluation of the partial 
rates. 

The first methods need an explicit construction of the final nuclear states in the context of 
a nuclear model (shell model, RPA, Fermi gas model etc.) and, therefore, they are, in general, 
more tedious. This difficulty is avoided in the approximations of the second category (see below). 



5.2.1 State-by-state calculation of the /z — e conversion rates 

With this method, by constructing explicitly the final nuclear states | /) in the context of a 
nuclear model, the branching ratio Rn e - is obtained by "summing the partial rates" for all 
possible excited states. By using the multipole expansion of the (//~,e~) operators Eqs. ( p9|) 
and (|30|) , the total rate matrix elements M 2 oi are evaluated by 

M 2 ot =E(2*+i)/ s 2 [£ fe) 2 Ei</«cii^ Mj n^>r+ 

S fexc I, J 

(S) 2 El^H^ ,S)J il^>| 2 ] (55) 

i,j 

(s=0 for the vector operator and s=l for the axial vector one). The first term in the brackets of 
Eq. (]55f) contains the contribution coming from all the excited states | f exc ) of the final nucleus 
(incoherent rate) and the second term contains the contributions coming from the ground state 
to ground state (coherent rate). 

State-by- state RPA calculations 

In actual shell model calculations it is quite hard to construct the final states explicitly in 
medium and heavy nuclei. In such cases one can employ the RPA approximation. In the context 
of the quasi-particle RPA the final states entering the partial rate matrix elements are obtained 
by acting on the vacuum | 0) with the phonon creation operator 

Q+{fJM) = Y\x { J^ T) A + {a, JM) - F a (/ ' J ' r) A{a,lM)] (56) 



(angular momentum coupled representation) i.e. | /) = Q + | 0). The quantities X and Y in Eq. 
fl5B| ) are the forward and backward scattering amplitudes. The index a, runs over proton (r = 1) 
or neutron (r = —1) two particle configurations coupled to J. 

The transition matrix element giving the partial rate I\-_^ from an initial state + to an 
excited state | /) takes the form 

(/ 1 1 f^ J 1 1 + > = £ Wi [ X^ugv}? + Y«Mvg>U$ } (57) 

The probability amplitudes V and U for the single particle states to be occupied and unoccupied, 
respectively, are determined from the BCS equations and the X and Y matrices are provided by 
solving the QRPA equations. The quantities W£ = W^ contain the reduced matrix elements 
of the operators T of Eqs. fl2"9|) and ( |30"D between the single particle proton or neutron states j± 
and ji as 

Wf 2jl = {j 2 \\f j \\j 1 )/(2J+1) (58) 



5.2.2 Effective computation of the individual /i — ► e partial rates 

The contribution of each individual state in the total fi — e conversion rate is effectively taking into 
account in the methods: i) nuclear matter mapped into finite nuclei and ii) sum rule approach 
(closure approximation). For a short description of these methods we devote the next subsections. 

1. Nuclear matter mapped into finite nuclei via a local density approximation 

This method, uses the exact muon wave function and gives the incoherent rate in terms of 



the Lindhard function 28 



where ni(q') and n2(q) the integral (0 or 1) occupation probabilities of proton-particle proton-hole 
(or neutron-particle neutron-hole), respectively. In the fi — e conversion this function corresponds 
to ph excitations of p-p or n-n type. 

This method resembles to the relativistic Fermi gas model and uses the local density approx- 
imation for the evaluation of the incoherent (/z - , e~) conversion rate, as 

r^-A - e~A*) = -2/^|$ M (x)| 2 /-^-n, 2?/; 



M n J (2tt) 3 l 2Ei 
x [SS|T| 2 (/i"p -^ e'p)ImU PtP (p lJL -p e ) 
+ SS|T| 2 (/i"n -^ e~n)ImU ntri (Pn ~ Pe)] (60) 

T represents the transition matrix, i.e. the amplitude associated with the elementary processes 
/i~p — ► e~p and jjTn —>■ e~n, and m^, Ei are the masses, energies of the particles involved in these 
reactions. The two terms in the brackets of Eq. (|60D result because of the charge conserving 
character of the (/i~,e~) operator. 



One can separate the nuclear dependence from the dependence on the elementary sector in 
Eq. (|B0D by factorizing outside the integral of eq. flBTf ) an average value of the quantity SS|T| 2 . 
For the ratio R of the incoherent (/i~, e~) rate divided by the total muon capture rate, we get 

R = roF^) = -T- ^ — 2 — J G( ' 4 ' z) (61) 

(the second step holds for the non-photonic case, involving the box diagrams of Fig. 1(c)). All 



the nuclear information is contained in the function G(A, Z) defined in Ref. | 2q| . Notice that 
R differs from the branching ratio R^e-, because R does not include the coherent part of the 
(/i~, e~) process. The calculation of R in various nuclei is based upon G(A, Z), since the \T\ 2 at 
the elementary level is in most models independent on the nuclear parameters. 

The branching ratio R^ e - in this method is obtained by adding in eq. ( |6T] ) the quantity 

R' = Vcoh ^ - 6 P (62) 

calculated independently. To calculate R' one can use the experimental data for the nuclear form 



factors (electron scattering, pionic atoms data etc. [|45| ) and the muon capture data [ftq] . 

2. Sum rule approach 

In the closure approximation, the contribution of each final state | /) to the total rate is 
also approximately taken into account. One assumes a mean excitation energy of the nucleus 
E = (Ef) —E gs , uses completeness relation for the final states | /), i.e. J2f \f)(f\ — 1> an d avoids 
the construction of each final state explicitly. 

The method proceeds by defining the operator taken as a tensor product from the single- 
particle operators T of Eq. ( p9|) or fl30l) . For a + initial (ground) state the relevant tensor 
product is 

= J2[T k xf k ']° Q (63) 



k,k' 

The corresponding total rate matrix elements are then written as 

ML = (^) 2 [ft(i I O v | i) + 3/i(i | 6 A | »)] (64) 

where | k | is the average momentum transfer. The operators Oy (vector) and Oa (axial vector) 
contain both one-body and two-body pieces. Consequently, in the closure approximation one has 
to evaluate the matrix elements of the operators Oy, Oa in the ground state \i) described by a 
specific nuclear model. 

The mean excitation energy E of the nucleus, involved in the matrix elements is defined 

as §4\ 

Y. f q 2 f\U\T J \gs)\ 

where the numerator represents the energy weighted sum rule and the denominator the non- 
energy weighted sum rule. 



Obviously, the mean excitation energy E should be evaluated by constructing explicitly all the 
possible excited states | /) in the context of a nuclear model. But in the closure approximation 
one wishes to avoid such calculations. Fortunately, E does not strongly depend on (A, Z) and 
its calculation for some properly chosen representative nuclei enables the use of closure for every 
nucleus throughout the periodic table. In Ref. p9| , the QRPA method has been used for the 



determination of the mean excitation energy for some nuclear isotopes including 48 Ti and 208 Pb 
which are of current experimental interest. 



6 Results and Discussion 



6.1 Nuclear structure dependence of the fi — > e conversion rates 

The nuclear physics calculations done for the (fi~,e~) conversion mainly refer to the most in- 
teresting quantities of the process, i.e. the branching ratio R^ e - and the ratio r\ of Eq. (|32j). 
Several nuclear methods for the most important mechanisms leading to this process have been 
employed. A general feature of the results obtained for various types of diagrams is the fact 
that qualitatively they are very similar (see as an example Fig. 3 below). The quantitative 
differences are mainly related to the presence or not of neutron contributions (they exist only 
for non-photonic diagrams) and the different values of the parameters /3,/3' and /3" of Eqs. (|13|), 
(TJ) and (|T5D which are the only elementary physics parameters entering the nuclear transition 



matrix elements. 

For the coherent rate and the total rate obtained using a sum rule approach, only the ground 
state of the nucleus need be constructed. The first such extensive calculations for closed (sub)shell 
nuclei were done with a determinantal ground state wave function (Slater determinant). These 
calculations provided a good description of the coherent rate p7| . The sum-rule results, however, 



have shown that the total rates are sensitive to the mean excitation energy required in the closure 



approximation. Furthermore, the local density approximation [^ gives reliable results for the 
coherent rate mainly due to the exact computation of the muon-nucleus overlap integral. The 
main advantage of this method in calculating the incoherent rate (by integrating over a continuum 
of excited states of a local Fermi sea), is the fact that it takes into account the contributions 
coming from the continuum, which are not explicitly included in RPA and shell model results. 

Below we discuss the results obtained in the framework of the quasi-particle RPA and discuss 
their comparison with those of shell model and local density approximation. The coherent and 
incoherent rate are discussed separately. 



6.1.1 The coherent process 

As we have mentioned in sect. 5, for the coherent rate the proton, F z (q 2 ), and neutron, F N (q 2 ), 
nuclear form factors are required. The results obtained by using quasiparticle RPA and shell 
model are shown in Fig. 3. They refer to the following two cases: (i) By neglecting the muon 
binding energy e& in Eq. ( |3~TD the elastic momentum transfer is the same for all nuclei, i.e. q ~ 
m 



/-< 



0.535/m^ 1 (results denoted QRPA(i)). (ii) By taking into account e& in Eq. fl3"T|), the elastic 
momentum transfer is q ~ m^ — e^ and varies from q « 0.529 fm~ l , for 48 Tz where e& ~ 1.3MeV, 
to q « 0.482/m™ 1 , for 208 Pb where e b « 10.5MeV (results indicated as QRPA(ii)). From the 



variation of the coherent nuclear matrix elements M^ oh with respect to A and Z illustrated in Fig. 
3(a) (photonic mechanism), and Fig. 3(b) (non-photonic mechanism), we see that by taking into 
consideration the muon binding energy e&, QRPA(ii), all matrix elements increase continuously 
up to the Pb region where they become about a factor of two larger than the approximate results 
of case (i). This implies that the coherent rate becomes larger for heavy nuclei, Pb region, which 
makes such nuclei attractive from an experimental point of view || ^TJ provided, of course, that 
they satisfy other additional criteria, as the minimization of the reaction background etc. ||. |5[ . 
The (fi~,e~) conversion electrons of a given target are expected to show a pronounced peak 
around E e = m M — ej, (for lead region E e « 95MeV). 



The above results have been obtained using the factorization approximation Eq. (54|). This 
approximation only slightly affects the ratio rj of the coherent rate to the total rate. However, the 
branching ratio R^ for heavy nuclei is appreciably affected by using exact muon wave function. 
This gives results which for heavy nuclei are 30 — 40% larger than the approximate ones. 

We should also note that, it is important to estimate the effect of the ground state RPA 
correlations on the (/i~,e~) matrix elements. This can be done by using a correlated quasi- 
particle RPA vacuum instead of the uncorrected one. In Ref. p9| such calculations have been 
performed using the correlated vacuum defined by K| E7 



| 0) = N e^ + \ 0) (66) 

where the operator S + contains the correlation matrix and Nq is an appropriate normalization 
factor. In this case, the coherent rate matrix elements take the form 

(0 | f | 0) = Nl (0 | f | 0) (67) 

which means that the correlated matrix elements are obtained by rescaling the uncorrelated 
ones (a similar expression holds for the total rate matrix elements). It was found that the 
matrix elements for 48 Ti obtained this way |29j are reduced by ~ 30% and that the ground state 



correlations tend to decrease the strengths of all (/x , e ) conversion channels. 



6.1.2 Incoherent process 

The results for the incoherent rate obtained with "state-by-state" QRPA calculations of all the 
excited states for a series of nuclei throughout the periodic table are shown in Table 2 for the 
photonic diagrams, and Fig. 4 for the non-photonic ones. Positive and negative parity states up 
to 6~, 6 + are included. For the photonic mechanism only the vector component, Sy gives non- 
zero contribution {Mf nc = Sy). For the non-photonic mechanism we have non zero contributions 
from both the vector and axial vector components, Sy and Sa, (M? nc = Sy + 3Sa)- 

It can be seen that, the main contribution to the incoherent rate comes from the 1~ mul- 
tipolarity. For every multipole only the low-lying excited states give significant contributions. 
High-lying excited states contribute negligibly. Furthermore, the incoherent matrix elements do 
not show clear A and Z dependence. Their magnitude depends on the spectrum of the individual 
nuclear isotope. 

In obtaining the results of Table 2 and Fig. 4 for 1~ states, the spurious center of mass 
contaminations, which arise from the use of empirical single particle energies and the truncation 
of the model space, have been eliminated by explicitly calculating the purely spurious state 



\S > and removing its admixtures from the incoherent and total rates. The calculation of the 
overlaps (l~,m\S) (where m counts the 1~ contaminated excited states) have shown that mostly 
the spuriousity lies in the lowest 1~ state being ~ 80% for 48 Ti and 208 Pb (for other nuclei 
see Ref. []53]). We showed that in a good approximation we can neglect the first 1~ state, but 
then we have to renormalize the strengths of the residual interaction (Bonn potential). From 
the renormalization we found that the incoherent \x~ — > e~ conversion rate is not sensitive to 
the strength parameters of the two-body interaction, which increases the reliability of the RPA 
results. For all nuclei studied, the spurious center of mass contribution is less than 30% of the 
incoherent matrix elements, i.e., ~ 2% of the total (fi~,e~) conversion rate. Furthermore, even 
after removal of the spurious components the 1~ multipolarity appears to be the most important 
incoherent channel. 

The explicit construction of all final nuclear states enables the calculation of the mean ex- 
citation energy E needed in the closure approximation. Though E is defined in analogy with 
the ordinary muon capture reaction by Eq. fl65|), the values of the "mean excitation energies" 
obtained in the (/i~, e~) are completely different than those evaluated in (jf , v^) [29| because in 



the last process the coherent channel doesn't exist. For the (f/,~, e~) process the coherent channel 
appears only in the denominator of Eq. QB5| ) and, since this dominates the total \i — e conversion 
rate, the resulting mean excitation energy E in this process is much smaller than that for the 
(/i - ,^) reaction (see Table 3). 



6.2 Comparison of coherent and incoherent processes 

As we have seen in sect. 4, a useful quantity for the (fi~,e~) conversion is the fraction rj of 
the coherent matrix elements M^ oh divided by the total one M t 2 oi . In earlier calculations rj was 
estimated to be a decreasing function of A with a value of r] ps 83 % in Cu region. By using, 
however, local density approximation [[2]J and QRPA |3l| we find that the coherent channel 
dominates throughout the periodic table. 



The values of 77 obtained with QRPA are a bit larger than the more accurate results obtained 
with a local density approximation [fjl| (LDA) in the medium and heavy nuclei region. The reason 
is the following: The LDA, for the part of non-coherent contributions, uses neither closure nor 
explicit summation over all final states. It makes a summation over the continuum of excited 
states (even though they contribute negligibly they are infinite in number) in a local Fermi sea. 
The incoherent rate obtained this way, is tied to the number of states one can excite; the larger the 
better; the muon mass provides the energy for the excitation of such states. Consequently, LDA 
takes into account the contribution of the continuum spectrum and calculates it more accurately 
for heavier nuclei. Such contributions are absent from the RPA results for every nucleus. 



6.3 Limits of the elementary sector part of R^-- Calculation of 7(-A, Z) 
The task of studying the nuclear physics aspects of the exotic (/z~,e~) conversion process is to 



evaluate the function 7 (A, Z) of Eqs. fl34j ), (pTj), ([450 an d fl47|). In general, the nuclear matrix 
elements entering j(A,Z) depend on the final nuclear state populated during process (3). The 
values for j(A, Z) calculated from the coherent QRPA matrix elements are shown in Table 4 



and compared with the results of Refs. p7| |28 |. The variation of j(A, Z) through the periodic 



table |27| exhibits a strong dependence on the neutron excess (N — Z) which mainly reflects the 



dependence on (A, Z) of the total muon capture rate. We should note that, the muon binding 
energy e& of Eq. Q5T|), which affects significantly the nuclear form factors for heavy nuclei like 
208 Pb, affects also the factor E e p e /m 2 in the definition of 7 Eq. fl34|). This factor takes into 
account the phase space in the transition matrix elements. By ignoring e&, this factor becomes 
unity. For i8 Ti this factor, by including e&, is equal to 0.97 and its neglect is a good approximation. 
For 208 Pb, however, this factor is equal to 0.81 which means that, in heavy nuclei the dependence 
on €b of the phase space cannot be ignored. 

By putting the experimental limits |5|, |S|] of the branching ratio R^- on Eq. fl55| ) we determine 
upper bounds on the parameter p for photonic, Z-exchange and W-box diagrams both in conven- 
tional extensions of the standard model as well as SUSY theories (for the currently interesting 
nuclei 48 Ti and 208 Pb see Table 4). We note that the values of p for 48 Ti are slightly improved 
over those of Table 2 of Ref. ||29|| , but those of 208 Pb are appreciably smaller than those of sect. 
5.1 of Ref. p7|. This big difference is due to the fact that, in the present work we have used 



the experimental limit of Ref. [g]. This experiment at PSI improved on the previous limit [21 
of 208 Pb, which was used in Ref. [J27| , by an order of magnitude. For both nuclei we found that, 
the limits obtained from various nuclear models do not significantly differ from each other. We 
should stress, however, that, the different limits for 208 Pb in the shell model results of Table 4, 
are due to the neglect of e& in Eq. ([31]) when calculating the nuclear form factors as mentioned 
above; its consideration gives similar results to those of LDA and QRPA for both mechanisms. 
This implies that all nuclear models studied here give about same values for p. 

One can use the limits of p to parametrize the muon number violating quantities entering 
Eqs. fl38D , (f40|) , (|46p and (|48|) and also to compare them directly to the value given from gauge 



models. As an example, we quote the value of p = 8.2 x 10 18 obtained in the supersymmetric 
model of Ref. |27] discussed above. This prediction of p is considerably smaller compared to the 



values listed in Table 4, which were extracted from experiment. 



7 Summary and Conclusions 

The main conclusions stemming out of the results discussed above can be summarized as follows. 
The coherent mode dominates throughout the periodic table indicated from the fact that rj > 
90%. 

The quasi-particle RPA and Local Density Approximation results for the (p~ , e~) rate do not 
show maximum around A « 60 as had previously estimated. This rate keeps increasing up to 
heaviest nuclei, i.e. region of 208 Pb, presently used at PSI as target. In evaluating the nuclear 
matrix elements the muon binding energy should not be ignored especially for heavy nuclei. 

The great part of the incoherent rate comes from the low-lying excitations. The 1~ state 
gives the maximum incoherent contribution even after removing the spurious admixtures. 

In contrast to the ordinary muon capture, the mean excitation energy, E, of the nucleus in 
pT — > e~ is very small, E « lMeV. 

The dependence of R^ e - on nuclear physics in coherent production, expressed by the function 
7 (A, Z) shows a large variation. Thus, for example, in the case of photonic mechanism, 7 lies in 
the region 1.6 < 7^ < 26.0 from lighter to heavier nuclei. 



All nuclear models used give about the same values for 7 and quantity p which contains the 
elementary sector parameters of the branching ratio R^ e -. The quantity p is useful to fix the 
lepton flavour violating parameters and test the various gauge models. 

The observation of any muon number (and in general any lepton and/or flavour number) 
non-conserving process, will reveal "new physics" beyond the standard model. The predicted 
branching ratios, however, are much smaller than experimental limits. The most optimistic 
results are obtained in the context of supersymmetry 1.2 x 1CT 18 < R^ e - < 2.4 x 1CT 16 

Although the predicted branching ratios for \i~ — ► e~ process is much smaller than present- 
day experimental sensitivity, this should not discourage the relevant experiments, since it is a 
common belief that we do not as yet have a complete theory to adequately describe such exotic 
processes. 

T.S.K. is grateful to the Organizing Committee of the Conference for financial support and the 
hospitality at Dubna. 
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Table 1. Best upper limits for the branching ratios of some muon number violating processes: 
Elementary particle decay modes and exotic neutrinoless muon decays in a nucleus. 



Process 


Best upper limit 


Reaction 


Ref. 


Elementary muon number violatin 


tg processes 




/i^e7 

\x — •> ee + e~ (// — > 3e) 

r ->//7 

r — > //e + e~, 

t — > e[i + fi~ 

T — > flfi + fl~ 

(/i+e-) <-* (/i-e+) 


R n < 4.9 x 10-n 

R 3e < 1.0 x 10-i 2 

R T < 4.2 x 10- 6 

R M m < 6.5 x 10~ 7 




[36] 
[37] 
[38] 

[39] 


Neutrinoless rare muon decays in 


a Nucleus 




H-(A,Z)-*e~{A,Z)* 

/2-(A,Z) ^e + (A,Z-2) 


R^- < 7.0 x IO-1 3 
R^ < 4.9 x 10-n 
R^+ < 5.5 x 10-12 


\x~Ti — > e~Ti 
irPb^e-Pb 
ji~Ti — > e + Ca 


[5] 
[6] 
[22] 



TABLE 2. Incoherent \x — e conversion matrix elements (Mf nc ) for the photonic mechanism. 
Only the vector component (Sy) of the /i — e conversion operator (see Eq. ( p7|) ) contributes. 



,r 


22 J 6 


™Ni 


HGe 


48 Gu 


162 Yh 
70 r ° 


IfPb 


0+ 


1.946 


1.160 


2.552 


2.088 


4.305 


2.512 


2+ 


0.242 


0.738 


1.396 


2.669 


6.384 


2.342 


4+ 


0.004 


0.005 


0.015 


0.021 


0.063 


0.056 


6+ 


6 10- 6 


6 10- 6 


1 10- 5 


6 10- 5 


2 10- 4 


3 10" 4 


1- 


3.711 


4.215 


5.066 


5.282 


4.824 


4.533 


3- 


0.037 


0.081 


0.152 


0.249 


0.542 


0.476 


5- 


2 10~ 4 


2 10" 4 


5 10- 4 


0.001 


0.005 


0.005 


M 2 

inc 


5.940 


6.199 


9.181 


10.309 


16.123 


9.924 



Table 3. Mean excitation energies for some representative nuclei in the reaction: /i^ + 
(A, Z) — > e~ + (A, Z)* . In contrast to the ordinary muon capture reaction, where E w 20MeV, 
in the ji — e conversion E is very smaller. 



Nucleus 


48 T* 


m Ni 


72 Ge 


U2 Cd 


W2y b 


208 pfo 


E (MeV) 


0.84 


0.44 


0.50 


0.40 


0.44 


0.28 



Table 4. The new limits on the elementary sector part of the exotic p — e conversion 
branching ratio extracted by using Eq. (|33| ) and the recent experimental data for the nuclear 
targets 208 Pb, 48 Tz, H [(J. The nuclear part of the branching ratio, described by the function 
j (A, Z) (see sect. 4), is also shown. 







48 Tl 


20£ 


] Pb 


Method 


Mechanism 


i{Az) 


p (xl0~ 13 


l(A,Z) 


p xl0~ 12 


QRPA 


Photon exch. 
W-boson exch. 
SUSY s-lept. 
SUSY Z-exch. 


9.42 

25.31 

25.61 

110.57 


< 4.6 

< 1-7 

< 1-7 

< 0.4 


17.33 

49.22 

49.50 

236.19 


< 2.7 

< 0.9 

< 0.9 

< 0.2 


LDA 


Photonic 
W-boson exch. 


9.99 
26.60 


< 4.3 

< 1.6 


17.84 
55.53 


< 2.6 

< 0.8 


SM 


Photonic 
W-boson exch. 


9.74 
26.50 


< 4.4 

< 1.6 


10.42 
27.43 


< 4.4 

< 1-7 



FIGURE CAPTIONS 



FIGURE 1. Typical diagrams entering the neutrinoless (fi~,e~) conversion: photonic 1(a), Z- 
exchange l(a),(b) and W-boson exchange 1(c), for the specific mechanism involving intermediate 
neutrinos. 



FIGURE 2. SUSY diagrams leading to the (/i~,e~) conversion: photonic 2(a), Z-exchange 
2(a),(b) and box diagrams 2(c), in a supersymmetric model with charged s-lepton and neutralino 
mixing are shown. Note that, the Z-exchange in 2(b) as well as in Fig. 1(b), comes out of 
electrically neutral particles (photon exchange does not occur in these diagrams). These Z- 
exchange diagrams may be important (those of Figs. 1(a) and 2(a) are suppressed by m^/m|). 

FIGURE 3. Variation of the coherent (fj,~,e~) conversion matrix elements M^ oh for specific 
mass A and charge Z (see text) for the photonic mechanism (a) and the non-photonic mechanism 
(b). In QRPA(i) the muon binding energy e& was neglected, but it was included in QRPA(ii). 
We see that e& strongly affects the matrix elements for heavy nuclei. For comparison the results 
of Ref. |27j (shell model results) are also shown. For photonic and non-photonic diagrams the 
coherent rate increases up to Pb region where it starts to decrease. 



FIGURE 4. Incoherent non-photonic (/i~,e~) conversion rate for the most interesting nuclei 
from an experimental point of view; 208 Pb is currently used as a target in the SINDRUM II 
experiment at PSI and A8 Ti has provided the most stringent limit on the flavour violation. The 
bars show the partial contribution of each multipolarity of the /i — > e conversion operator (black 
bars for the vector part (Sy) an d empty bars for the axial vector part (<SU))- The 1~ contribution 
shown is obtained after removing the spurious contaminations but still it is the most important 
one for this process. 
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